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Abstract 


Aerodynamic drag on moving vehicles consists the major cost element of 
transportations. There have been recent developments, where suppression of wake 
by the use of edge treatment have promised dramatic savings in fuel costs. 

In present work a ’full vortex cloud’ scheme has been used as a tool for nu- 
merical investigation of unsteady, incompressible flows past some two dimensional 
bodies. Four geometries have been considered, first geometry is circular cylinder 
for validation of the code, second one is a square cylinder to be used as a model 
for trucks, third is a similar square cylinder with edge modifications, and forth is 
a rotating circular cylinder. The drag obtained for square cylinder are compared 
with the available results, and it is found that the results are in good agreement 
with the available results It is also found that the the drag of square cylinder can 
be reduced to a good extent by using a drag suppressor. We have also tried to 
develop the flow past rotating circular cylinder. 
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Chapter 1 


Introduction 


1.1 Literature survey 

The origin of vortex methods can be traced back to 1930s in the works of 
Rosenhead(1931) who approximated the motion of a two dimensional vortex sheet 
by following the movement in time of a system of point vortices, using their loca- 
tions as quadrature points. Performing his calculations by hand, he was not able 
to conduct simulations for extended times. 

Westwater(1935) was first to apply the approach of Rosenhead to the vortex 
sheet roll-up behind a wing. With the advent of high speed computers, a number of 
researchers e.g.Takami(1964), Moore(1974), Clements and Maull(1973) extended 
this approach to the vortex sheet roll-up behind a wing. Some researchers Kuwa- 
hara and Takami(1973) and Chorin and Bernard(1973) introduced vortices with 
small but finite area instead of using point vortices. 

In 1973 Moore presented a numerical study for the roll-up of finite vortex sheet. 
In his work he introduced tip vortex to represent the tightly rolled portion of vortex 
sheet, the chaotic motion which was a feature of some of the earlier studies, was 
eliminated and details of the outer portion of the spiral were calculated. 

In 1973 Chorin presented a scheme to simulate viscous effects. In his scheme 
vortices move under mutually induced velocity of field modified by a diffusion 
component due to viscosity. In his method of random walk particles undergo a 
Brownian-like motion to simulate the effects of viscosity. 

In 1979 Baker used a numerical method to calculate the time evolution of vortex 



sheet by adapting ’’Cloud m Cell” technique. Krasny in 1987 studied numerically 
the two vortex sheet evolution problems arising in the aerodynamics(e.g. aircraft 
wake) . 

Other researchers worked for the development of vortex methods include Sarp- 
kaya(1989), Smith and Stansby(1988), Saffman(1981), Milizzano and Saffman(1977) 
and Leonard(1980) etc. 


1.2 Present Work 

The present work is in essence an extension of the work started in 1993 by 
Parag Kumar continued by Verma(1994) and Manish Kumar Jain(1996). Parag 
Kumar studied the motion of the vortices shed from sharp corners by a method of 
conformal transformation. Jain and Verma worked on operator splitting method to 
solve the Navier-Stokes equation and developed a full- vortex cloud method. They 
solved the flow past circular and square cylinder and provided basic technique for 
flow past multiple objects. 

The present work is focused on developing the technique further to study com- 
plex shapes. 



Chapter 2 


THE FLOW PAST SQUARE 
CYLINDER 


2.1 Introduction 

Vortex methods were originally conceived as tools to model the evolution of un- 
steady, incompressible, high Reynolds number flows of engineering interest. Vortex 
methods simulate flows of this type by discretizing only the vorticity carrying re- 
gions and tracking the computational elements in a Lagrangian frame. 

The ’’full- vortex cloud theory” attempts to develop the real flow past two- 
dimensional bodies of arbitrary shape. Vorticity is created over the whole body 
surface, i.e. body is covered with discretized vortex panels. The strength of which 
is determined using Martensen method which is based on satisfying the zero tan- 
gential velocity condition by postulating a vortex sheet at the body surface. There- 
after this vortex sheet is diffused in the flow field and is convected downstream as 
a cloud of discrete vortices during a sequence of small time steps, ultimately shed- 
ding naturally from sharp edges or boundary layer natural separation points. The 
ultimate aim is nothing less than a full simulation of real fluid flow with minimum 
human intervention. 



2.2 The full- vortex cloud approach 


In this method, the traditional discrete vortex methods are modified so as to 
account for viscous diffusion as well. This permits calculations for finite Reynolds 
numbers as opposed to traditional methods which are valid only for infinite Reynolds 
number. The surface of two-dimensional body is covered by discretized vortex pan- 
els. These panels then diffuse into the flow region with every time step as well as 
convect according to the local flow velocity. 

The complex diffusion convection problem is handled by an operator-splitting 
strategy first proposed by Chorin(1973, 1978) and later developed extensively by 
Smith and Stansby(1987), Lewis and Porthouse(1983) and Verma(1994). 

Thus the Navier-Stokes equation 

d S + q.Vq = ~ + vV 2 q (2.1) 

is first converted into the two dimensional vorticity transport equation 

^ + q.Vco = iA7 2 w (2.2) 


and then split into two equations 


• convection 


du 

~di 


+ q.Vu — 0 


(2.3) 


and 

doo o / . 

• diffusion — = vV 2 uj (2.4) 

The above two equations are solved one after the other. The process starts with 
the calculation of potential flow about the body with the vortex sheet replacing 
the body surface. This sheet is then discretized into a number of elements. These 
elementary vortices first diffuse according to equation 2.4 and, then, are convected 
according to equation 2.3. The convection velocity is determined by the influence 
of all the vorticity elements, the free-stream of the new vortex sheet required to 
satisfy tangency condition. After these two equations have been solved, the Navier- 
Stokes equation reduces to 

| (2.5) 

dt p 

which when solved will give pressure distribution. 


The numerical scheme is developed in this chapter. The flow diagram in Fig. 
2.1 shows the strategy for the full- vortex cloud analysis. 




Figure 2.1: Flow diagram for full-vortex cloud model 


















2.3 Martensen method 


Consider a discretized two-dimensional body which, is immersed in a uniform 
flow Wqo at an angle of attack measured with respect to the x-axis as shown 
m Fig 2 2. Let the length of n th panel be ds n and has vorticity sheet of strength 
t( s ) P er un it length Then from Biot-Savart law, the induced velocity at s m , dq mn 
due to a small vorticity j(s n )ds n located at s n on the body is 


d tmn = (2,6) 

where r mn is the distance between the elements m and n. 


dm 



Figure 2 2: Martensen method 

Resolving dq mn parallel to s m , we obtain the following tangential velocity com- 
ponent 

, 1 f (ijm ~ yn) c O s 0m ~ { x m ~ X n )sin(3 m ^ f fr) 

**“ = 2? I (*.-*.)’ + <!*.-*)» J 7(s " )<is “ (%1) 

Stating the Dirichlet boundary condition at s m 

— 2 T ^ k(,s m , s n ) , j{s n )ds n T (UooCOSpm + Wo sinfim) 0 (2-8) 

Coupling coefficient k(s m , s n ) is given by 

1 f ( Urn ~~ yn) C0S Pm ~ ( x m ~ £n) 5W /?Tn 1 

{Sm ’ Sn> ~ to \ (i„ -x n ) 2 + (y m - y n ) 2 I 


( 2 . 9 ) 



Coupling coefficient K(s m ,s n ) can be written as 


S n ) Stl ) (4 10 ) 

In special case when n = m, the coupling coefficient can be modified as follows, 

(2.11) 

where K mm can be evaluated using L— Hospital’s rule 


]_ f 

K(s m , s m ) = — - + K mm 


ry-' y J (i/m Vn^COsPrn (*£m *^n) SITlPn 

J. x fyyi />nn _ 11X11 


or 


2tT Sm->$n 1 ( “1“ ( Dm Un 

tv -' ^ ^-Pm 

mm ~ 47rr m ” ~4tT 

Lewis gave a simple relation for self induced coupling coefficient, 

K(s m j S m ) ~ 2 g^(An+l Pm— l) 

Equation 2.8 can be written m the following numerical form 

M 


n=l 


(2 12 ) 


(2 13) 


(2 14) 


m? ) n f{s n ) — (UooCOsPm + VooSVTlpm) (2.15) 


Martensen method in presence of cloud of vortices 


In presence of the cloud of vortices, the boundary integral equation will yield, 

-\li s m) + j K S m, S n )7( s n)d5 n + ( U^COsPm + V^sinfim) 


+ ^ ^ AT j (UmjCOsPm ”1“ VmjSlTlj3 m ) 0 

3=1 


(2.16) 


or m numerical form 


M 


K(s m , s n )q(s n ) = —(UooCOsfim + Voo sinpm) — ATjlUmjCOSfim + VmjSinPm ) 

n = 1 

(2.17) 

Where the unit velocities U m) , V mJ are given by 


Uma ~ 2 tt 


Vm Vj 


m3 


(2.18) 



(2.19) 


Where r 2 m] = (x m - x 3 f + (y m - Vj ) 2 

Helmholtz theorem states that the vorticity at any moment in the flow should 
be conserved. Thus 


M Z 

Y 7(Sn) A S n + Y^J~ Tare = 0 (2.20) 

n=l j=l 

Where Y arc is the sum of strength of all the vortices removed from the flow field. 
The above equation (2 11) must be added to equation (2.11) to ensure vorticity 
conservation. Thus 

M 

y>(*ffl, S n ) + As n ) 7 (s n ) = -(UocCOsPm + VooSinPm) 

n = 1 


z 

^ ' ^Y 3 (U m jCOS f3 m T V m jSlTl/3 m + 1) + T clrc 

. 7=1 


(2.21) 


2.4 Input data and data preparation for body 
geometry 

The body surface is discretized into small elements, the numbering of the ele- 
ments is done in clockwise direction starting from the leading edge. The last point 
on the body coincides with the first point. The profile slope and profile pivotal 
points are calculated at the middle of each element. 

A s n = yj (z„ + i - x n ) 2 + (y„+ 1 - y n ) 2 (2.22) 

profile slopes and pivotal points which are located at the centre of each element, 
are then given by 

cosf}„ = ~ . ? - 4 (2.23) 

/Xs n 

Sinp n = ^ n+ A X (2.24) 

x n = ^ ( x n+l + x n) (2.25) 

Vn - + w) (2.26) 



Input data points 



Figure 2.3: Surface element geometry 


2.5 Calculation of coupling coefficients 

The calculation of coupling coefficient is straight forward except for some com- 
plications in special circumstances, such as for thin bodies. For thin body shapes, 
the coupling coefficients are unable to represent correctly the mutual interference of 
opposite elements correctly. For thick body shapes on the other hand this problem 
does not arise and the coupling coefficient matrix is diagonally dominant. 

While in case of thin bodies, the error level of back diagonal coefficient of the 
matrix, namely K(s m ,SM+i-m), may dominate leading to erroneous results. To 
eliminate the error the following techniques are used. 

1. Use of sub-elements 

A standard Martensen numerical model implies the replacement of the line 
vorticity elements by a concentrated vortex of strength 7 (s n )ds n at each pivotal 
point x.e. at the middle of each element. 

The representation of the true element vorticity influence can be improved when 
K(s m , s n ) is calculated as the average value of coupling coefficients for elements. 

Consider a case in which an element s n is broken-up into N sub-elements each 
of length ds n /N. The coupling coefficient of the m th element is then given by 

K(s m , s n ) = — / k{sm)S n )d$ n (2,27) 

JQ 



In numerical foim, 


K^Sm-, &n) 


2 ttN 


N 

E 

Z =1 


{'Urn Videos ^Yn {x^n Xi'jSlTljSn 

(x m ~ X t ) 2 + ( y m ~ Vi ) 2 


(2.28) 


Where the location of i th sub-element is given by 

x, = I. + - i(l + 1 V)~j As «^ s i 3 ’‘ (2,29) 

y. = v« + (» - ;(i + ao) (2.30) 

It is found that two or three sub-elements give better results than using more sub- 
elements. 


2. Back-diagonal correction 

Back diagonal correction is used to enforce zero internal circulation From 
Kelvin’s theorem the net circulation Ar TO around the profile interior induced by 
surface vorticity element 7(s m ) As m should be zero. If this condition is forced upon 
the coupling coefficients, 

1 M 

K(s opp , s m ) — ~r ^ K(s n , s m )As n ( opp — M 1 tti) (2.31) 

n=l, n^opp 

The back diagonal matrix coefficients K(sM+i-m, s m ) are replaced by the value 
given by this equation, thereby, ensuring net circulation around the profile interior 
implied by the numerical model is made to zero. 


2.6 Shedding of vortices from the body 

Following potential flow analysis at each time step, the surface vorticity 7(s„) As n 
created at every boundary element is shed from the surface into the fluid as a new 
discrete vortex The vorticity is released at some offset from the body. The offset 
is taken as y/iiTKt/S. However at very high Reynolds number, the offset is very 
small. The usual approach is to choose a value for the offset related to element 
size(say 25% of the average surface element length). 



2.7 Diffusion of vortices 


As given by Chorin m 1973, the diffusion of vortices in a flow is modeled 
by a random walk method. The diffusion model is such that the vortices move 
randomly in ^-direction with flat probability, and walking in the r-direction with 
normal distribution 

Let P denotes the probability that a given vorticity element will lie somewhere 
within the circle of radius r, 


p = 1 - e ~ r2 l Avt 

(2.32) 

for the random walk process, we generate a random number Pi between 0 to.l with 
uniform distribution. The i t h vortex element will be found at radius such that 

-r? 

P z = 1 — e 

(2.33) 

the radial shift of vortex can be calculated from the above 
\Jkvt ln(l/Pi). We generate two random numbers P l and Q % . 
step At, the element i will undergo the following movement, 

equation as r % = 

Thus during time 

• angular shift A 8 t = 2nQ z 

(2.34) 

• radial shift Ar t = ^ kv At\n(l j Pf) 

(2.35) 

Thus after the increment of time step At the new location of the i th element will 

be 


x[ = x l + Ar l cosA8 l 

(2.36) 

y[ = yt + Ar,smA6 l 

(2.37) 


2.8 Convection of vortices 

After shedding and diffusion having been taken place, vortices undergo convec- 
tion under the influence of all other vortices in the flow plus vorticity bound to the 
body surface plus the external flow field. 

If the m th vortex is at ( x , y) and there are total Z vortices in the flow field, then 
the convection velocity of the vortex is the sum of individual velocities imparted 




Figure 2.4’ Random number geneiator 


by each vortex on this particular vortex, plus sum of velocities imparted by the 
bound vorticity of the body plus free stream velocity. In numerical form, 

V, — Uqo "4" Uyortex cloud d~ U bound ( 2 

V = V V ortex cloud d~ Abound (2.39) 

where = W^cos® and Voo = W^sma 

If we consider vortex in the cloud located at the position (xj . yj ), the velocity 
at m due to vortex n will then be given by Biot-Savart law. Summing up the 
influence of Z vortices as 


u , * - Y" — ( - (2.40) 

u vortex cloud / 2^ \ r 2 / 

(2.41) 

If there are M discrete panels on the body surface with distributed vorticity of the 
strength 7 (s n )ds n at the n th element, 

tt 7( 5 j)^ s j f v 

U bound / j 2^ l 7*2 

— 1 \ 



(2.42) 






(2.43) 


V b 


bound 


Y' 7( 3 j)^gj f 

27r l r 2 ) 

j=i v 7 


where r = \/{x — ^) 2 + (y — y^) 2 

Now after convection velocity of each vortex has been calculated, their new 
position after convection is ascertained by displacing them by amounts Ax = 
uAt and Ay = vAt and adding these displacements to their original location 
coordinates will give their new location coordinates^. e 


x =x + uAt (2.44) 

y = y + v At (2.45) 

After convection of all the vortices, if we now re-calculate the convection ve- 
locities (u ,v) at their new locations (x ,y), we can estimate the better location 
by using the central difference scheme. That is 


x = x + \-{u + u)At 

(2.46) 

v" = y + ^(v + v')At 

(2.47) 



(a) Forward difference method 


(b) Central difference method 


Figure 2.5: Estimates of streamline path 


2.9 Merging of vortices 

As the vortices move away from the body, their effect diminishes with increasing 
the distance from body, we know that CPU time required for the N particles in 
the flow is of the order 0{N 2 ). So reducing the number of vortices by merging 
them is the effective way to reduce their number. The total strength of the newly 
formed vortex after merging will be the algebraic sum of the strength of individual 


vortices. 



2.10 Reflection of vortices 


Vortices if drifted very close to the body, can produce serious potential flow 
errors So vortices are pushed back in the direction normal to the body to some 
fixed distance (say offset e) in order to prevent the possible errors in Martensen 
analysis or in convection 


2.11 Calculation of lift and drag coefficient and 
Strouhal number 


Since we are using operator splitting technique for the solution of Navier-Stokes 
equation, after diffusion vV 2 q and convection q.Vq, the remaining part of the 
equation is, 


Vp _ dq 
p dt 


( 2 . 48 ) 


At any point s n on the body surface, the velocity q = 7 (s n ) parallel to the surface. 
Then the pressure gradient along the surface at s n is given by, 


dp s 

ds 


~P 


d jjsn) 

dt 


which can be written in the following form, 

j(s n )ASr, 


A p n = -p- 


At 


AT n 
) At 


( 2 . 49 ) 


( 2 . 50 ) 


for pressure change over the surface element n during the discrete time step At. 
The integration of above equation will give pressure change at any element on the 
body surface. 

m 

Pm=Pl+J2 A Pn ( 2 - 51 ) 

n= 1 

The average pressure on anv element will then be, 


Pn 


Pn + Pn + 1 \ 

2 J 


( 2 . 52 ) 


After calculating average pressures on all the elements, the lift and drag coefficient 
can be calculated by the following equations, 


Ci = 




M 

^ ^ Pn C0S PnAs n 
n=l 


( 2 . 53 ) 



(2.54) 


C A = 


~ TT/2 / ^ ^ Pn Sm !^n-^ s r. 


kpKV 

where l is the characteristic length. 

To determine the Strouhal number St of vortex shedding, spectral analysis is first 
used to determine the dominant frequency of the pressure data. Then the dominant 
frequency is used to determine the St using the following equation 


(2.55) 


where d is the diameter of the cylinder, and lh T 0O is the free-stream velocity. 



2.12 Results and Discussions 


We have first run the code for the case of circular cylinder for the purpose of 
validation of the code. The results obtained for circular cylinder are then com- 
pared with the available results The Fig. 2.6A shows the available numerical 
results(j?e/.ll, pp. 422) and Fig. 2.6B shows the results obtained in the present 
study. The Fig. 2.6A predict the average lift and drag of circular cylinder as zero 
and 1.5 respectively at Re=20000. And from the present study, the lift and drag 
obtained are zero and 1.53 respectively at Re=20000. Thus the results obtained 
are in fair agreement with the available numerical results. After validation, the 
flow past a square cylinder at various Reynolds number is calculated. 

Flow past square cylinder 

To solve flow past square cylinder, the side of square cylinder is taken as 1.0, 
mainstream flow W, x is 1.0 and a oo = 0. Here we have discretized the body into 56 
elements. The time step is taken sufficiently small at 0 05 to maintain reversibility 
during convection and the code is run for a minimum of 650 steps(T mm = 32.5) 
which is sufficiently large to cover several oscillations of the vortex sheet. The 
lift and drag coefficients along with Strouhal number are calculated at various 
Reynolds number(20, 000, 30,000, 40,000, 50,000, 60,000, 1x10 s , 1.2 xlO 6 ). For 
eliminating numerical noise, we are averaging the lift and drag coefficients for 10 
steps. The flow diagram for different Reynolds numbers shown in Fig. 2.8 to Fig. 
2.12. The variation of lift and drag coefficients with respect to time are shown 
from Fig. 2.13 to Fig. 2.15. The Table 2.1 shows the average values of life and 
drag coefficients with Strouhal number. 

The average drag obtained from the present study for square cylinder is around 
2.18 at Re = 1 x 10 5 . The experimental results(Ref. 25) show the drag of square 
cylinder around 2.10. Thus the results obtained numerically are in good agreement 
with the available results. Fig. 2.7 shows the variation of drag with increasing 
Reynolds number. The drag although slightly higher than the experimental results, 
is in good agreement with the experimental results. 



Lift and drag 



Time 


Figure 2.6B:Vortex cloud solution obtained for flow past circular cylinder at Re = 
20000 for validation of the code 
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Re 

Average lift 

coefficients, Ci 

Average drag 

coefficients, Cd 

Strouhal No. 

St 

20000 

0.011405 

2.299080 

0.1953 

30000 

-0.054072 

2.154425 

0.2013 

40000 

0.137076 

2.113487 

0.2106 

50000 

-0 015132 

2 328437 

0.2216 

60000 

-0.011620 

2.177026 

0.1973 

1 x 10 5 

0.127729 

2.177774 

0 1946 

1 2 x 10 6 

-0.098907 

2.190457 

0.1920 


Table 2 T Lift and Drag for the case of square cylinder 
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Variation of drag of square cylinder with Reynolds number 



Reynolds number 


Figure 2.7: Comparison of drag of square cylinder with available experimental 
results(i?e/. 25) 





























Lift and drag Lift and drag 
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Variation of lift and drag of square cylinder at Re=30000 



Figure 2.13: Lift and drag variation at Re = 20000 & 30000 for square cylinder 
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Lift and drag Lift and drag 
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Vanation of lift and drag of square cylinder at Re=60000 



Figure 2.14: Lift and drag variation at Re = 40000 & 60000 for square cylinder 





Lift and drag Lift and drag 


Variation of lift and drag of square cylinder at Re=l*10 A 5 



Variation of lift and drag of square cylinder at Re=1.2*10 A 6 



Figure 2.15: Lift and drag variation at Re = lxlO 5 & 1.2 xlO 6 for square cylinder 





Chapter 3 


THE FLOW PAST SQUARE 
CYLINDER WITH 
ATTACHMENTS 

3.1 Introduction 

Flow separation on a body prevents complete recovery of pressure and hence re- 
sponsible for drag. Drag forces are large in case of bodies having sharp edges. So 
the idea here is to delay the onset of separation by using attachments with the 
square cylinder. Ritter[16] and Som[21] have proposed techniques, using simple 
attachments to reduce drag on open railway wagons and on truck bodies. Experi- 
mental investigations have shown significant reduction in drag. 

We propose here to study some geometries to see if the full- vortex cloud method 
can predict reduction in drag. 


3.2 Mathematical Formulation 

To solve the flow past this type of arrangements, we develop the flow past multiple 
bodies first. The difference in the flow past single bodies and the multiple bodies 
lies in the determination of coupling coefficient matrix. In latter case, the effect 
of a body over the other is also taken into account while calculating the coupling 
coefficient matrix. Rest of the treatment remains the same. Suppose there are P 



bodies in x — y plane immersed m uniform flow 11'^ at an angle of attack with 
respect to x— axis. The body q is discretized into M q number of elements with 
number of bodies from 1 to P. 

Thus the Boundary-Integral equation taking into account the effects of one body 
over the other can be written as, 

p Mq 

y(Sqn) U (Xi COs3p m I -c SlTl3p m (3 1) 

q= 1 n=l 

where also also p = 1, 2, P and m — 1,2, ..., M p 

The above equation 3 1 is the statement of the Dirichlet boundary condition at 
the element m of body p. The coupling coefficient i.e. the velocity induced on one 
body element due to the element of the other body is given by, 


K pq 


^Sq n ( ( Upm VqnjCOS Ppm {%pm 


2tt 


(3.2) 


(Zpm %qn)“ T (l/pm Vq-ri)' 

n) is the self induced coupling coefficient 


For special cases when (p = q and m 
as given earlier for single body. 

The equation 3.1 can be written in the following matrix form, 


^■11 A12 

A 2 1 A 2 2 


[ 7 (s)] = [rhs] 


(3 3) 


Here sub-matrices An and A 22 are the coupling coefficient matrices as calculated 
earlier for single body. And sub-matrices A i2 and i4 2 i contain the effect of one 
body over the another. Typically, sub-matrix A\ 2 contains all of the coefficient 
accounting for interference experienced by the first body due to the second body. 


Body -2 



Body -2 



Figure 3.1: Two models used in present investigation 





3.3 Results and discussions 


We have considered two geometries. In both of the geometries, the side of square 
cylinder is taken as 1 0, free stream velocity is 1.0, angle of attack, a M = 0, 
length of the flat portion of the attachment is 0.5, thickness and inner radius of the 
curved part of the attachment is 0.1 and 0.4 respectively. Both of the geometries 
are studied at different positions of attachment and for various Reynolds number. 


For first geometry 

The geometry is shown in Fig. 3.2. The body of square cylinder and the attach- 
ment are discretized into 56 surface elements and 52 surface elements respectively. 
The Table 3.1 shows the average values of lift and drag coefficients of the main 
body(square cylinder) and the attachment for various gaps between the square and 
the attachment at Reynolds number of 20000. Total lift and drag acting on the 
geometry is also shown in the Table 3.1. It is evident from the Table 3.1 as the gap 
is reduced the total drag acting on the geometry reduces and reaches to a mini- 
mum, and on further reduction of the gap, drag increases. Maximum reduction in 
drag is around 38% at Reynolds number of 20,000. 


Attachment (Body -2) 



Figure 3.2: First geometry 


The case is again run for various Reynolds numbers(at 50, 000, 1x10 s and 1.2x 
10 6 ) when the gap between the two bodies is optimum(i.e. when the drag is min- 
imum). The results are shown in the Table 3.2. The flow diagram for various 



Gap 

Average lift coefficients, Ci 

Average drag coef ficients, Cd 

Reduction 


Body — 1 

Body — 2 

Total 

Body — 1 

Body — 2 

Total 

in C d {%) 

0.20 

-1.078648 

0.135818 

-0.942830 

1.548154 

0.111806 

1.659960 

27.80 

0.17 

-1.194436 

0.173199 

-1.021237 

1.494855 

0.117757 

1.612612 

29.86 

0.16 

-0.644665 

0.192763 

-0.451902 

1.239422 

0.186909 

1.426331 

37.96 

0.15 

-1.427683 

0.186493 

-1 241190 

1.732633 

0.047211 

1.779844 

22.58 

0.14 

-1.116135 

0.160531 

-0.955604 

1.671501 

0.134620 

1.806121 

21.44 


Table 3.1: Lift and drag of geometry-1 at various positions of the attachment 


gaps are shown in Fig. 3.6 to Fig. 3.9. Fig. 3.10 shows the flow diagram for 
Re = 1.2 x 10 6 at optimum gap. The variation in lift and drag coefficient with 
time is shown Fig. 3.11 to Fig. 3.12. The Fig. 3.13 to 3.14 shows the varia- 
tion in lift and drag at optimum gap and Re=l x 10 5 and 1.2 x 10 6 respectively. 
The present study shows a reduction in drag using drag suppressor from 2.18 to 
lA28{reductionof 38%) at Reynolds number of 20000. 


Re 

Average lift coef ficients, Ci 

Average drag coefficients, Cd 

QHH 


Body — 1 

Body — 2 

Total 

Body — 1 

Body — 2 

Total 

■ 

Ifrfr • 


-0.644665 

0.192763 

-0.451902 

1.239422 

0.186909 

1.426331 

37.96 

50000 

-1.716633 

0.197500 

-1.519133 

1.471544 

0.070148 

1.541693 

33.79 

lx Id 5 

-1.072953 

0.159764 

-0.913189 

1.435434 

0.140764 

1.576198 

27.62 

1.2 x 10 6 

-1.103621 

0.207633 

-0.895989 

1.376179 

0.142850 

1.519029 

30.65 


Table 3.2: Lift and drag of geometry- 1 at optimum gap 


For second geometry 

The geometry under consideration is shown in Fig. 3.3. The body of square 
cylinder and attachment are discretized into 56 and 38 surface elements respec- 
tively. The Table 3.3 shows the lift and drag coefficients of individual bodies and 
the total lift and drag on the geometry at Reynolds number of 20,000 when the 

gap between the bodies vary from 0.14 to 0.20. It can also be seen from the Table 
3.3 that the drag acting on the geometry is minimum when the gap is 0.16 and 






















increases when the gap is either reduced or increased. 



Figure 3.3: Second geometry 

The case is again run for the configuration of minimum drag for various Reynolds 
numbers(at 50, 000, 1 x 10 5 and 1.2 x 10 6 ). The results are shown in the Table 3.4. 
Flow diagrams for various gaps are shown in Fig. 3.15 to Fig. 3.18. Fig. 3.19 
shows the flow diagram for Re — 1.2 xlO 6 and optimum gap. The variation in lift 
and drag with time is shown in Fig. 3.20 to Fi.g 3.23. Fig. 3.24 to Fig. 3.25 shows 
the variation in lift and drag of the geometry with time when the gap is kept at 
optimum and Re = 1 x 10 5 and Re = 1.2 x 10 6 respectively. 


Gap 

Average lift coefficients, Ci 

Average drag coef ficients, Cd 



Body — 1 

Body — 2 

Total 

Body — 1 

Body — 2 

Total 






1.723650 

0.071530 

1.795180 

21.90 




-1.009274 

1.490483 

0.024947 

1.515430 

34.08 

0.16 

1.610553 

0.268493 

- 1.342060 

1.344694 

0.014013 






-1.624638 




37.00 

0.14 

-1.660850 

0.257079 

-1.403771 

1.582787 

0.038369 

1.621156 

29.48 


Table 3.3: Lift and drag of geometry-2 at various positions of the attachment 


Fig. 3.4 shows the variation in drag of geometries when gap between the bodies 
is varied and average drag of square cylinder is also plotted. The plot shows that 














Re 

Average lift coef ficients, Ci 

Average drag coe f ficients. C d 

Reduction 


Body — 1 

Body — 2 

Total 

Body — 1 

Body — 2 

Total 

inC d {%) 

20000 

-1.610553 

0 268493 

-1.342060 

1.344694 

0.014013 

1.358707 

40.90 

50000 

-1 365549 

0 239934 

-1.125615 

1.421387 

0.016249 

1.437636 

38.26 

lx Id 5 

-1 383413 

0.230505 

-1 152908 

1.441183 

0.046595 

1.487778 

31.68 

1.2 x 10 6 

-1.565719 

0.276271 

-1.289448 

1.498898 

0.004638 

1.503535 

31.36 


Table 3.4: Lift and drag of geometry-2 at optimum gap 


the drag of cylinder with attachments is less than that of square cylinder alone 
The drag of both the geometries is minimum when the gap is 0.16. 

In Fig. 3 5A drag for square cylinder and attachments for both of the geometries 
are plotted for increasing Reynolds number. In the same plot, drag of square 
cylinder is also plotted. In Fig. 3 5B, total drag acting on the geometries are 
compared with the drag of square cylinder. It is clear form the plot that drag 
reduces to a good extent when using drag suppressor. 


Comparison of drag of model- 1 and model-2 at varying gaps with square 



Figure 3.4: Comparison of drag of geometries at various gaps between the two 


bodies 
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Comparison of drag of square cylinder with and without attachments 
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Figure 3.5: Comparison of drag of square cylinder with both of the geometries 




Figure 3.6: Flow diagram for geometry-1 at Jle = 20000 and gap = 0.20 between 
the two bodies 


















Figure 3.9: Flow diagram for geometry-1 at Re = 20000 and gap = 0.14 between 


the two bodies 








Figure 3.10: Flow diagram for geometry-1 at Re = 1.2 xlO 6 and optimum gap 


between the two bodies 









Lift and drag Lift and drag 


Variation of lift and drag of geometry- 1 at gap=0.20 



Variation of lift and drag of geometry-1 at gap=0.17 



Figure 3.11: Variation of lift and drag of geometry- 1 at varying gaps between the 
bodies at Reynolds number 20000 





Lift and drag Lift and drag 


Variation of lift and drag of geometry- 1 at gap=0.16 



Variation of lift and drag of geometry- 1 at gap=0.14 



Figure 3.12: Variation of lift and drag of geometry-1 at varying gaps between the 
two bodies at Reynolds number 20000 


















Figure 3.15. Flow diagram for geometry-2 at Re = 20000 and gap - 0.20 between 


the two bodies 


















Figure 3.19: Flow diagram for geometry-2 at Re = 1.2 xlO 6 at optimum gap 


between the two bodies 

















Variation of lift geometry-2 at gap=0.16 







Variation of lift geometry-2 at gap=0.14 



Variation of drag of geometry-2 at gap=0.14 



Figure 3.23: Variation of lift and drag of geometry-2 at Reynolds number 20000 
and gap= 0.14 between the two bodies 
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V anation of lift of geometry-2 at gap=0. 16 
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Variation of drag of geometry -2 at gap=0.16 



Figure 3.24: Variation of lift and drag of geometry-2 at optimum gap between the 
two bodies at Re = 1 x 10 5 








Chapter 4 


THE FLOW PAST ROTATING 
CIRCULAR CYLINDER 

4.1 Problem Formulation 

Flow past rotating circular cylinder can be modeled by introducing a bound vor- 
tex F at the centre of cylinder. The bound vortex T can be related to the bound 
vorticity by, 

r = j) 7 (s)ds 

or in numerical form, 

M 

F = 

71=1 

and Martensen’s boundary integral equation for plane 

M 

n)7(Sn) = -(UooCOSPm + V^Sinfim) (4.3) 

n= 1 

Adding equations (4.2) and (4.3), the modified form of Martensen equations is, 

M 

^ ] ( K(s m , S n ) + As n ) 7 (s n ) = ( Uqo cosfim + VaQSinPm) + r (4.4) 


(4.1) 


(4.2) 

two dimensional flow is, 


In presence of cloud of vortices the equation 4.4 is modified as 

M 

^(^( 5 m) s n) + ^Sn)j(Sn) ~ —{U^COSj^m + \ ^ Siri{3 m ) +F 

n=l 

Z 

“i“F circ - Y2 ^j(Um]COsP m + \ mjSmdm + 1) (4 5) 

3=1 

where U mj andV m;j are the same as given earlier for square cylinder. 

The bound vortex T can be related to the rotation ft by F = 27 r rfi ft. Thus 
the only change for solving the flow past rotating circular cylinder is change in 
Martensen equation and the rest of the analysis(i.e. shedding, diffusion, convection 
etc.) remains unchanged. 


4.2 Results and discussions 

The diameter of cylinder is taken as 1.0, mainstream flow is 1.0, and rotation(Q) 
of cylinder is varied. The cylinder is discretized into 32 surface elements The flow 
is calculated for 720 time steps, each time step is 0.05. We also solved the flow for 
infinite Reynolds number(Potential flow) . The result for infinite Reynolds number 
and for vortex cloud solution are given in Table 4.1 The results obtained for lift 
and drag are not correct. Flow diagrams for different rotation speeds are given in 
Fig. 4.2 to Fig 4.6. Variation in lift and drag with time is shown in Fig. 4.7 to 
Fig. 4.11. The Fig. 4.1 shows the variation in lift and drag of the rotating cylinder 
with increasing velocity ratio(rO/C/oo). 



Rotation 

At Infinite Reynolds 

1 ’ ortex cloud analysis. 

n 

number (Potential flow ) 

at Re = 

:lxl0 5 


Ci 

c d 

Ci 

c d 

0.1 

0.313943 

0.000000 

-0.218555 

0.643247 

0.2 

0.627886 

0 000000 

-0.612226 

0.421329 

0.4 

1.255773 

0 000000 

-0.654632 

0.715022 

0.6 

1 883659 

0 000000 

-0.759954 

0.575640 

0.8 

2.511545 

0.000000 

-0.901063 

0.395768 

1.0 

3.139432 

0 000000 

-0.103302 

0.799423 

1 2 

3.767318 

0.000000 

-0.688595 

0.548724 

1.4 

4.395204 

0.000000 

-0.589837 

0.647782 

1 6 

5.023091 

0.000000 

-0.506253 

0.398093 

1.8 

5 650977 

0.000000 

-0.446093 

0.612781 

2.0 

6 278863 

0.000000 

-0.105586 

0 748433 


Table 4.1: Lift and drag of rotating circular cylinder at varying rotation speeds 
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Lift and drag of rotating cylinder at increasing velocity ratios 



Figure 4.1: Lift and drag of rotating cylinder at increasing velocity ratio, W( 


















































Chapter 5 


CONCLUSIONS 


In first part of the present work we developed a procedure for the flow past 
square cylinder in flow field. The lift and drag coefficients are calculated by ob- 
tained pressure distribution on the surface of the cylinder, and are compared with 
the available experimental results. The drag and lift coefficients obtained are in 
fairly good agreement with the available experimental results. 

The second part of the thesis deals with the extension of flow past square 
cylinder to flow past square cylinder with drag suppressors. We first developed 
a general procedure to solve the flow past two-dimensional multiple bodies. Two 
geometries of drag suppressors are dealt in with and it is found that drag of square 
cylinder can be reduced to a large extent, as has been reported in the literature. 

In the last part of thesis, we tried to solve the flow past rotating circular 
cylinder. In the available time, this part could not be taken to a satisfactory 
conclusion and the results obtained are not correct. 
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